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This article addresses a modification of local time for stochastic processes, to be 
referred to as natural local time. It is prompted by theoretical developments arising 
in mathematical treatments of recent experiments and observations of phenomena in 
the geophysical and biological sciences pertaining to dispersion in the presence of an 
interface of discontinuity in dispersion coefficients. The results illustrate new ways in 
which to use the theory of stochastic processes to infer macro scale parameters and 
^ behavior from micro scale observations in particular heterogeneous environments. 
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> 1 Introduction and Motivating Example 

5^ The purpose of this article is to present some new results that call attention to the presence 

and effects of particular types of heterogeneity observed to occur in diverse geophysical and 
ecological/biological spatial environments. A few examples of interest to the authors are 
provided to highlight the occurrence of interfacial discontinuities, but one may expect that 
readers will easily be able to conceive of many others. Apart from the motivating examples, 
theoretical results are provided to illustrate the interfacial effects in ways that may prove 
useful in analyzing and interpreting more complex data sets. 
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The most basic empirical considerations involve various temporal measurements, e.g., 
breakthrough times, residence or occupation times, recurrence times and, as will be ex- 
plained, a fundamental notion of local time. 

Mathematical local time is a quantity with a long history in the theory of stochastic 
processes. It concerns a mathematical measure of the amount of time a stochastic process 
spends locally about a point, which we will refer to as mathematical local time in distinction 
with a modification introduced here to be referred to as natural local time^ Starting with 
the celebrated theorem of Trotter [IH] establishing the continuity of mathematical local time 
for standard Brownian motion, the general problem of determining necessary and sufficient 
conditions for continuity of local time for stochastic processes is an old and well-studied 
problem at the foundations of probability theory; see [I7j for recent progress and insight into 
the technical depth of the problem for a large class of Markov processes. 

While the present paper has little to offer to the mathematical foundations of the subject, 
observations from field and laboratory experiments are provided together with some theo- 
retical results that point to new directions in this area from the perspective of applications 
and modeling of certain dispersive phenomena. At the most fundamental level, the issue 
involves the units of measurement. In particular, the units of mathematical local time in the 
context of dispersion of particle concentrations are typically those of spatial length. This is 
not unnatural in the mathematical context owing to the locally linear relationship between 
spatial variance and time when the dispersion coefficient is sufficiently smooth. However, as 
will be seen, this breaks down in the presence of an interface of discontinuity in the diffusion 
coefficient. 

Effort is made to make this article accessible to a diverse audience by including some 
basic mathematical definitions and background theory. A complete and systematic treat- 
ment of essentially all of the underlying mathematical concepts can be found in [H]. We 
conclude this introductory section with an empirical example that serves to motivate and 
illustrate the role of a key mathematical construct, namely that of skew Brownian motion. 
This is followed by a section with additional empirical illustrations from the ecological and 
biological sciences. The emphasis of the paper is on mathematical theory and modeling. The 
main results demonstrate (i) the effects of general point interface transmission probabilities 
on breakthrough times and (ii) on residence times; (iii) symmetry (via martingale) rela- 
tionships between interfacial transmission probabilities, dispersion coefficients, and skewness 
parameters; (iv) a special role for continuity of flux (conservation of mass) verses continuity 
of derivatives at an interface; (v) a special role for (natural) local time continuity in the 

^We reserve the use of the term "physical local time" for a more specialized case, but physical and/or 
biological modeling considerations underlie the development of the nomenclature introduced in this article. 
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determination of the transmission parameter at the interface. 

Example A (Dispersion in Porous Media). The topic addressed in this paper origi- 
nally initiated as a result of questions resulting from recent laboratory experiments designed 
to empirically test and understand advection-dispersion in the presence of sharp interfaces, 
e.g., experiments by [23], [20], [Z]- Such laboratory experiments have been rather sophisti- 
cated in the use of layers of sands and/or glass beads of different granularities and modern 
measurement technology (see Figure [T]). 
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Figure 1: Saturated Fine/Coarse Grain Sand Separated by an Interface Experimental 
Assembly, Oregon State University 

From an engineering science point of view, basic research in this area is aimed to improve 
understanding of the role of sharp discontinuities in the hydrologic parameters appearing 
in equations used for predictions of the spread of contaminants in saturated porous media. 
There is a huge literature in the sciences and engineering pertaining to the general problem. 
For a perspective on the mathematical foundations see [TU], [H],[S], [IH], [33], [S]- The former 
considers heterogeneity from a deterministic framework, while the latter involve randomness 
in the medium. p3j contains a survey of various approaches, together with a more exhaustive 
set of references to the geophysical and hydrologic literature. 

From a general mathematical point of view an interface is defined by a hypersurface 
across which the dispersion coefficient is discontinuous. As is well-known for the case of 
dilute suspensions in a homogeneous medium (e.g., water), perhaps flowing at a rate v, the 
particle motion is that of a Brownian motion with a constant diffusion coefficient D > 
and drift v. When applied to colloidal suspensions, this is the famous classic theoretical 
result of Albert Einstein [16] experimentally confirmed by Jean Perrin [H]; see [22] for a 
more contemporary experiment. The essential physics of the problem, however, is relatively 
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unchanged when one considers dissolved chemical species. 

The basic phenomena of interest to us here is captured by the following: 
Question. Suppose that a dilute solute is injected at a point L units to the left of an 
interface at the origin and retrieved at a point L units to the right of the interface. Let D~ 
denote the (constant) dispersion coefficient to the left of the origin and that to the right, 
with say < (see Figure [T]). Conversely, suppose the solute is injected at a point L 
units to the right of the interface and retrieved at a point L units to the left. In which of 
these two symmetric arrangements will the immersed solute most rapidly breakthrough at the 
opposite end ? 
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Figure 2: Interfacial Schematic 

In the experimental set-up, the coordinate in the direction of flow is one-dimensional flow 
across an interface. As will be explained, a localized point interface results in a skewness 
effect that explains much of the empirically observed results suggested above; see [37], [38] . 

13], m m, m- 

Let us briefly recall the notion of skew Brownian motion introduced by Ito and McKean 
in ^2]J to identify a class of stochastic processes corresponding to Feller's classiflcation of 
one-dimensional diffusions. The simplest deflnition of skew Brownian motion B^"'^ for a 
given parameter, referred to as the transmission probability, a G [0, 1] is as follows: Let \B\ 
denote reflecting Brownian motion starting at zero, and enumerate the excursion intervals 
away from zero by Ji , J2 , . . . . 

This is possible because Brownian motion has continuous paths and, therefore, its set of 
zeros is a closed set; the set's compliment is open and, hence, a countable disjoint union of 
open intervals. Let Ai,A2,... be an i.i.d. sequence of Bernoulli ±1 coin tossing random 
variables, independent of B, with P{An = 1) = a. Then, B^"^ is obtained by changing the 
signs of the excursion over the intervals J„ whenever to y4„ = — 1, for n = 1, 2, . . . (see Figure 
[T|. That is, 

00 

(1.1) i?(")(t) = ^A„U(t)|i?(t)|, t>0. 

n=l 

In particular, the path t — )■ B^°'\t),t > is continuous with probability one since the 
Brownian paths are continuous. 
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Figure 3: Skew Brownian Motion Construction 

For positive parameters D^, D~ , consider a piecewise constant dispersion coefficient with 
interface at a: = given by 

D{x) = D-l(_oo,o)(a;) + D+l[o,oo){x). xeR, 

Here Ij denotes the indicator function of an interval J defined by lj{x) = l,x & J, lj{x) = 
0,x ^ J. In response to the question about breakthrough times raised above, the following 
theorem provides an answer as a consequence of [37], [3]. For convenience of mathematical 
notation take f = for here. 

Remark 1.1. Just as standard Brownian motion may be obtained as a limit of rescaled 
simple symmetric random walks, e.g., fTl^ - an alternative description of skew Brownian 
motion with transmission probability a may be obtained as a limit in distribution of of the 
rescaled skew random walk defined by ±1 displacement probabilities, 1/2 each, at nonzero 
lattice points, and having probabilities a, 1 — a for transitions from to +1 and to —1, 
respectively, see ll9j . 

Theorem 1.2. Let D^, be arbitrary positive numbers, with < . Define Y^°'*\t) = 
(j{b\" ^),t > 0, where B'^^*^ is skew Brownian motion with transmission parameter a* = 
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Vz)+ Wj*- ' """^ ^^^^ ^ A/D+xl[o,oo)(a;) + V^D~" a;l(_oo,o](a;), G I^- Let Ty = mi{t > 
= y}. Then, 

(a) For smooth initial data cq, c{t,y) = EyCQ{Y^°'\t)), t>0, solves 



dt 2dy dy ' dy dy 

(h) For y>Q, P.yiTy > t) < ^Py{T^y > t) < Py{T^y >t), t> 0. 

Observe, for example, that by integrating the complementary distribution functions in 
(b), one obtains that the mean breakthrough time in fine to coarse media is smaller than 
that of the mean breakthrough time a coarse to fine media. Related phenomena and results 
on dispersion in this context, including the case v ^ 0, are also given in |37], [38], [2], [3], 
[39], [M]- The proof of part (b) of the theorem relies on a transformation to elastic skew 
Brownian motion to eliminate the drift. In addition, a recently obtained formula for the first 
passage time distribution for skew Brownian motion is given in ||5j . Considerations of both 
local time and elastic standard Brownian motion appear in [15] to cite another biological 
context for their significance. 

The identification of the stochastic particle motions in the presence of an interface for 
simulation can have utility in problems involving the computation of particle concentration at 
a single spatial location, e.g., for so-called resident breakthrough, since other pde numerical 
schemes generally involve computation of the entire concentration curve. Some illustrative 
results pertaining to Monte-Carlo simulations of skew diffusions are described in [^ and 
references therein. In addition, the identification of skew diffusion answers the basic physics 
question of finding the particle motion that Jean Perrin would have reported had there been 
an interface ! 

From the macro-scale quantity of particle concentrations, the determination of a* in 



Theorem L2 may be viewed as the result of mass conservation. The interface condition 
is simply continuity of flux at the interface. However, as will be illustrated by examples 
in the next section, it is not necessary that continuity of flux always be obeyed in the 
presence of interfaces. Moreover, even the modeling may be at the scale of a single particle, 
e.g., biological dispersion of individual animals, in which 'particle concentrations' are not 
relevant, e.g., see [22], [2S], [12] for such examples. A stochastic particle model plays an 
essential role in such situations. 
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Figure 4: Fishing vessels off tfie coast of Argentina, http://subdude- 

"site.com/WebPics/WebPicsEarthLights/south-america-hilights-429x600.jpg courtesy of 
Ricardo Matano, COAS, Oregon State University 



2 Related Examples from Other Scientific Fields 

As illustrated by the examples in this section, the role of interfacial phenomena is of much 
broader interest than suggested by advection-dispersion experiments. However the specific 
nature of the interface can vary, depending on the specific phenomena. We briefly describe 
below three distinct classes of examples of phenomena from the biological/ecological sci- 
ences in which such interfaces naturally occur, together with compelling questions having 
substantial biological implications. 

Example B (Coastal Upwelling and Fisheries). Up-wellings, the movement of deep 
nutrient rich waters to the sun-lit ocean surface, occur in roughly one percent of the ocean 
but are responsible for nearly fifty-percent of the worlds fishing industry. The up-welling 
along the Malvinas current that occurs off of the coast of Argentina (see Figure |4]) is unusual 
in that it is the result of a very sharp break in the shelf, rather than being driven by winds. 
The highlighted points in the figure represent a flotilla of fishing boats concentrated along 
the shelf where the up-welling occurs. The equation for the free surface elevation r/ as a 
function of spatial variables (x, y) is derived from principles of geostrophic balance and takes 
the form 

, dr] r f dh\~^ d^T] 

dy f \dx J dx"^ ' 

where r > and / < in the southern hemisphere, and h{x) is the depth of the ocean at a 
distance x from the shore. In particular, the sharp break in the shelf makes h'{x) a piecewise 
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constant function with positive values ^ H~ . The location of the interface coincides with 
the distance to the shelf-break. In particular, if the spatial variable ?/ > is viewed as a 
'time' parameter, then this is a skew- diffusion equation, but the physics imply continuity of 
the derivatives drj/dx at the interface rather than 'flux' (see [23 cind references therein). 

Example C (Fender's Blue Butterfly). The Fenders Blue is an endangered species of 
butterfly found in the pacific northwestern United States. The primary habitat patch is 
Kinkaid's Lupin flower (Figure |2]). 



Fender's Blue 



Kincaid's Lupin 



Patch Distribution 





Figure 5: Fender's Blue Butterfly Habitat - Kincaid's Lupin Patches - U.S. Fish 
and Wildlife Service, Wikipedia http://en.wikipedia.org/wiki/Fender ?o27s_blue_butterfly, 
http : // www. fws. gov / oregonfwo /Species /PrairieSpecies / gallery, asp 

Quoting IHj, "Given past research on the Fender's blue, and the potential to investigate 
response to patch boundaries, we ask two central questions. First, how do organisms respond 
to habitat edges? Second, what are the implications of this behavior for residence times?" 
Sufficiently long residence (occupation) times in Lupin patches are required for pollination, 
eggs, larvae and ultimate sustainability of the population. Empirical evidence points to a 
skewness in random walk models for butterfly movement at the path boundaries. The deter- 
mination of proper interface conditions is primarily a statistical problem in this application, 
however we will see that certain theoretical qualitative analysis may be possible for setting 
ranges on interface transmission parameters. There is a rapidly growing literature on the 
statistical estimation of parameters for diffusions, [2S], |13], [I]- However much (though not 
all) of this literature is motivated by applications to dispersive models in finance for which 
the coefficients are presumed smooth and the data is high frequency. 

Example D (Sustainability on a River Network). The movement of larvae in a river 
system is often modeled by advective-dispersion equations in which the rates are determined 
by hydrologic / geomorphologic relationships in the form of the so-called Norton laws. In gen- 
eral river networks (Figure |2]) are modeled as directed binary tree graphs and each junction 
may be viewed as an interface. There are also special relations known from geomorphology 
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and hydrology that can be apphed to narrow the class of graphs observed in natural river 
basins; e.g., see |12], [M], [H], and references therein. 




Figure 6: Digitally Generated River Network Using HidroSIG Java; courtesy of Jorge 
Ramirez 

Conservation of mass leads to continuity of flux of larvae across each stream junction as 
the appropriate interface condition. Problems on sustainability in this context are generally 
formulated in terms of network size and characteristics relative to the production of larvae 
sufficient to prevent permanent downstream removal at low population sizes; see [lO] for 
recent results in the case of a river network, and [3D] and [TH] for related mathematical 
considerations. 



3 Natural Dispersion and Natural Occupation Time 

The following theorem provides a useful summary in one-dimension of the interplay between 
diffusion coefficients and broader classes of possible interfacial conditions. To set the stage 
we begin with Brownian motion X = {X{t) = \/DB(t) + vt : t > 0} with constant diffusion 
coefficient D > and constant drift v; i.e., B = {B{t) : t > 0} is standard Brownian motion 
starting at zero, with unit diffusion coefficient and zero drift. The stochastic process X has 
a number of properties characteristic of diffusion. Namely, the random path t — )■ X{t) is 
continuous with probability one, and 

(3.1) E{X{t)\X{u),u< s) = X{s) + v{t- s), 0<s<t. 

This makes X a continuous semimartingale. In the case f = 0, X is in fact said to be a 
(continuous) martingale. In general a continuous semimartingale is a process with continuous 
paths that differs from a (local) martingale by a (unique) continuous (adapted) process 
having finite total variation (in the sense of functions of bounded variation). Such structure 
is a natural consequence of the interfaces discussed in this article. This will be made clear 



at equation (3.12) below 
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The Brownian motion X is also a Markov process with (homogeneous) transition prob- 
abihties 

(3.2) p{t,x,y) = -j^^e-'-^^^, t>0, 

satisfying the advection-dispersion equation 

dp ^ d D dp dvp 
^ ■ ^ dt dy^2 dy' dy ' 

with continuous derivatives of all orders everywhere. Since the displacements X{t) — X{s) 
are independent for disjoint intervals (s,t), it follows that the variance of X{t) is linear in t 
and, in particular, E(X(t) — ft)^ = Dt,t > 0. The definition of mathematical local time at 
a, denoted if (a), for X can be expressed as dif{a) is the amount of time X spends in the 
infinitesimal neighborhood (a, a ± da) prior to time t. More precisely, 

(3.4) £f(a)=lim^/ l(„+,)(X(s))rfs. 

Note the presence of the diffusion coefficient D. Thus, in the context of dispersion of particles 
in a fluid, for example, [D] = L'^/T and [e] = [X] = L yields units [if] = j^x^xT = L. 

The standard extension of the mathematical definition of local time for a continuous 
semimartingale X exploits the quadratic variation < X >t,t > 0, of the process. If X is 
a square-integrable martingale then < X >j is defined by the property that the process 
X^(t)— < X >t,t > 0, is a martingale; see [H]. So, in the case of the Brownian motion 
one may check that < X >t= Dt,t > 0. In the case of a continuous semimartingale, the 
quadratic variation is quadratic variation of its martingale component. 

Mathematical local time at a can be defined as an increasing continuous stochastic process 
C(t),t > such that 



(3.5) 



|X(t) -a| = |X(0) -a| + /" sgn{X{s) -a)dX{s) +if{a), t > 0, 

Jo 



where sgn{x) = x/\x\,x ^ 0, and s(y'n(0) = —1 (by convention). For purposes of calculation 
it is often convenient to consider (right and left) one-sided versions defined by 

(3.6) £f '+(a) = lim - / (X(s))c? < X >„ t > 0, 

40 e Jo 

(3.7) £f'-(a) =lim- / l(„](X(s))rf<X>„ t > 0. 

40 e Jo 
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Then 

(3.8) £f (a) = * ' ^ ^ t > 0. 
The utihty of these quantities rests in the following two formulae: 

Ito-Tanaka Formula: If X is a continuous semimartingale with local time £^(a), a G M, t > 
0, and / is the difference of two convex functions then 

(3.9) /(X(t)) = /(X(0)) + f r_{X{s))dX{s) + if'^ia)f"ida), 

Jo ^ Jr 

where f"{da) is a positive measure corresponding to the second derivative of / in the sense 
of distributions. 

Occupation Time Formula: For a non-negative Borel function g and t > 0, one has a.e. 

that 

(3.10) ! g{X{s))d < X >s= [ g{a)de^{a). 

Jo Jr 

Theorem 3.1. Let D^,D^ be arbitrary positive numbers and let < a, X < 1. Define 
y("^(t) = a{B\"'*),t > 0, where B^"'^ is skew Brownian motion with transmission parameter 
a and a{x) = \/D+xl[o,oo){x) + Vd^ xl(_oo,o](a^)5 x G M. Then 

M{t) = f{Y^''\t)) -^-I^D{Y^-\u))nYiu))du, t > 0, 

zs a martmgale for all f eVx = {f e C2(M\{0}) n C(M) : A/'(0+) = (1 - A)f (Q-)} if and 
only if 

a = aiX) = — ; ; . 

Av^+(1-A)v^ 

Proof. In the case of skew Brownian motion one has the following relationships between 
mathematical local time at zero and its one-sided variants, e.g., see [31] 



(3.11) if''''\0) = 2aif"\0), £f'"''"(0) = 2(l-a)£f"'(0), t>0. 
Moreover, B^°'^ is the unique strong solution to the stochastic differential equation, see 

(3.12) dS(")(t) = ^^d£f"^'^{0) + dB{t). 

2a 



2a-l rtB(") 



In particular, considering the integrated version, B is the martingale component and 
is the finite variation component in the view of B^"'^ as a semimartingale. 
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It is also straightforward to relate mathematical local times of Y^"'"' and B^°'^ through the 
one-sided formulae as: 

£f "^+(0) = v^£f' "''+(0), eJ'"''-io) = v^if'"''-{o) 

Applying the Ito-Tanaka formula to the positive and negative parts of Y^°'^ together with 



(3.12), one has 



1 



and 



[B(«)(t)<0] 



dB{t) - 



{2a - 1) , 



2a ' ' ' 2 

Thus, since F*^"^ is the difference of its positive and negative parts and noting the cancellation, 
one has 



^^^^ + VTF'^]dir'^io). 

2 2a ' ^ ^ ' 



(3.13) dY^°'\t) = \ D{B^^)(t))dB{t) + 



For a difference of convex functions / G Pa? 



f"{da) = f"{a)da + (/'(0+) - f'{0-)6o{da). 



With these preliminaries, again use the Ito-Tanaka theorem together with (3.11) and (3.13) 
to get 

/(y(")(t)) = /(F(")(0))+ f nY^-\s))dY^-\s) + \ I C'^\a)nda) 

Jo ^ Jr 



/(F(-)(0)) + ^^ f'_iY^^\s))^D{Yi-){s))dB{s) 



+ 



av^ - (1 - a)v^- 
2a 



+ I I il^"^^^{a)f"{a)da 



/i(F(°)(s)X^"^'+(0) 



(3.14) 



+ 2 / ^r"''na)(/'(0+)-r(0-)5o(cia). 



According to the mathematical occupation time formula relating mathematical occupation 
time and mathematical local time, and and noting that the quadratic variation is given by 
d < >,= L>(r(")(s))rfs, it also follows that 



r(aX^"^'+(a) = I f D{Y('^\s))nY(^\s))ds. 
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The asserted result now follows by subtracting this term from each side of the equation (3.14) 



and noting the cancellation of local time terms leaving a stochastic integral with respect to 
Brownian motion when and only when a = a (A), i.e., 



aVD+ - (1 - a)VD- 
2a 



2a 



/'(0-)£f ^'+(0) + -v^ff •+(o)(/'(o+) - /'(o-)) 
i^v/D^/r)]£r'+(0)=0 



if and only if ^^"^^ = That this is enough follows from standard theory, e.g.. Theorem 



2.4 in 



□ 



Some implications for the examples will be described below in terms of the stochastic 
particle evolution, however one may also note that one has the following consequence at the 
scale of Kolmogorov's backward equation. 

Corollary 3.2. Let D^,D^ be arbitrary positive numbers and let < a, X < 1. Then for 
Co G Vx, the unique solution to 

- = -D(x) — 
dt 2^ 'dx^- 

is given by 



Oc dc 
^hm c(t,a;) =Co(x), A — (t,0+) = (1 - A) — (t, 0"), t > 0, 



It is illuminating to consider Theorem 3.1 in the context of the examples. 



Example A. Theorem 3.1 provides a generalization of the results obtained in [37] and [3] 
for the case of advection-dispersion problems across an interface described in Example A. 
One may check that 



(3.15) 



A 



D- 



D+ + D- 



a 



ZD- 



follow from Theorem 3.1 for this application. This coincides with the results of [37] and [3] 
obtained by other methods. 

The following definition is made with reference to both the diffusion coefficient and the 
interface parameter in the context of this and the other examples. 



Definition 3.3. With the choice of a = a{X) given by Theorem 3.1, we refer to the pro- 



cess y'^"'^'^)) as the natural diffusion corresponding to the dispersion coefficients D^,D and 
interface parameter A. 
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Additional Nomenclature: We sometimes refer to the natural diffusion corresponding to 
A = T^+2r,- ; ^* — hw^Ht^ physical diffusion. We refer to the diffusion in the case 



A = l/2,a^ = -1==^-^= = 1 — a* as the Stroock-Varadhan diffusion since it is the solution 



D++D- ' " y^+VD 

1/2 = ^ = 

to the particular martingale problem originating with these authors; see 

Example B. Observe that in the application to the coastal up- welling problem one obtains 



(3.16) A = 1/2, a(l/2) = = — 

VD+ + VD- 

The interface parameter provides continuity of the derivative, rather than of the flux, at 
the interface. The natural diffusion for this example may be checked to coincide with the 
Stroock-Varadhan diffusion in this case. Note that the answer to the first passage time 
problem will be exactly opposite to that obtained for the advection-dispersion experiments 
of Example A under this model. 

The following modification of the usual notion of mathematical occupation time, where 
the integration is with respect to quadratic variation and in units of squared-length, provides 
a quantity in units of time that we refer to as natural occupation time. |^ 

Definition 3.4. Let X be a continuous semimartingale. The natural occupation time of G 
by time t, is defined by f{G,t) = l(X(s) G G)ds,t > 0, for an arbitrary Borel subset G 
ofR. 

The following result illustrates another way in which the issue raised in Example B 
relating interfacial conditions to residence times is indeed a sensitive problem. The proof 
exploits the property of skew Brownian motion that for any t > 0, 

(3.17) P(fi(") > 0) = a. 

This is easily checked from definition and, intuitively, reflects the property that the excursion 
interval Jn{t) of \B\ containing t must result in a [An(t) = +1] coin flip, an event with 
probability a. 

Theorem 3.5. Let y denote the natural diffusion for the dispersion coefficients D^, D~ 
and interface parameter A. Denote natural occupation time processes by 

r+{t) = r 1[f("W)(s) > 0]ds, t > 0. 
Jo 

^ While we emphasize "natural" choices from the point of view of modeling (and units), there are very 
sound and important reasons for the standard mathematical definitions. In particular, no suggestion to 
change the mathematical definition is intended. Indeed, as the proof of Theorem |3.1| demonstrates, the 
notion of mathematical local time and occupation time and their relationship is extremely powerful in 
singling out the special value of a (A) for given interface parameter A and dispersion coefficients D^. 
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Similarly let (t) = t — r^(t), t > 0. Then, 
Ef+(t) > Ef^{t) Vt > 



A > 



with equality when A 



/D+ 



/D++VD- 



Proof. Using the definition of natural diffusion for tlie parameters i5^,A and tlie 

above noted property (3.17) of skew Brownian motion, one lias 

ft 



Ef+(t) := E [ l[r("(^))(s) > 0]ds 
Jo 

= e[ 1[v^5("(^))(s)) > 0]cis 
Jo 

= [ P(5("(^))(s) > 0)ds 
Jo 



= ta{\). 
Similarly Ef^(t) = t(l - a(A)). Thus 

Ef+(t) _ ^ «(A) 



ET-{t) 'l-a(A)' '(l-A)v^ 



The assertion now follows. 



□ 



Example A vs Example C. As noted previously, in the Example C pertaining to insect 
movement, the determination of an interface condition is largely a statistical issue as there is 
no scientific rationale to apply mass conservation principles, or smooth Fickean flux laws. In 
fact, it is interesting to observe that under the mass conservation one arrives at the interface 
parameter 

(3.18) X = D+/{D+ + D-), 

for which the residence time is longer on average in the region with the faster dispersion rate! 
While this is to be anticipated for physical experiments of dispersion in porous media of the 



type described in Example A, Theorem 3^ shows that in fact the conservative interface con- 
dition (defined by this choice of A) would not be appropriate for models of animal movement 
for which the faster dispersion occurs in more hostile environments; e.g., see [25], [12], [32] 
for relevant considerations. There is indeed something to be learned from data as to what 
exactly might apply, but such theoretical insights can provide a useful guide, and help to 
prevent mistaken assumptions when transferring more well-developed physical principles to 
biological/ecological phenomena. 
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4 Continuity of Natural Local Time 



We now discuss another issue pertaining to the definition of mathematical local time and 
the basis for the suggested modification to natural local time. 

In the often cited article jlT], it was first observed that skew-Brownian motion provides 
a non-trivial example of a continuous semimartingale on the interval (—00, 00) having a 
discontinuous mathematical local time. Let us now examine this situation in the context of 
natural diffusions. For the purposes of this discussion take the drift v = and consider a 
constant diffusion coefficient = D~ = D 

Since the quadratic variation of skew Brownian motion coincides with that of Brownian 
motion, one has < \/DB^°'^ >=< \/DB > . Among this class of natural diffusions, one 
may ask what distinguishes the particular diffusion X = \/DB ? Of course, the answer is 
that X is determined by a = A = 1/2. If one views this choice in the context of the flux 
in particle concentration, then it provides continuity of flux. On the other hand, in view 
of the respective theorems of Trotter [46j and Walsh [57], it is also the unique choice of a 
from among all skew Brownian motions to make local time continuous. The latter may be 
viewed as a stochastic particle determination of the physical diffusion model, among natural 
diffusions, for constant diffusion coefficient D. 

The next theorem, a version of which was originally conceived in [4j , extends this to the 
more general framework of the present paper, in particular to include the case 7^ D~ . 
However, it requires the following modification of the definition of mathematical local time, 
referred to here as natural local time. 

Definition 4.1. Let X be a continuous semimartingale. The natural local time at a if = 
- — 2"^^ — of X is defined by 

(4.1) £f'+(a) =lim- / l[„,,+,)(X(s))ds, £f (a) = lim - / l(,_,,,](X(s))ds 

40 e Jq 40 e Jo 

providing the indicated limits exist almost surely. 

The units of natural local time are then [if] = ^, appropriate to a measurement of 
(occupation) time in the vicinity of a spatial location a. While the purpose here is not to 
explore the generality for which natural local time exists among all continuous semimartin- 
gales, according to the following theorem it does exist for natural skew diffusion. Moreover, 
continuity has a special significance. 

Theorem 4.2. Let Y^"'^^^^ be the natural skew diffusion with parameters D'^,X. Then the 
natural local time ofY^"'^'^^^ at is continuous if and only if X = j^+'^jj- , i-^-, if o,nd only if 
a(A) = a* and thus F^"*) is the physical diffusion. 



16 



Proof. In view of (3.11) observe that 



(0) 



a 



1 — a' 



Vt > 0. 



Moreover, 



lim - 

40 e 



1[0 < < e]ds 

1 \/r)+ /"* 

= lim^^ / 1[0 < < 



'D 
1 



Similarly, (0) 



'"''-(0). Thus, 



-(0). 



'D 



a(A) 



z]ds 



Now simply observe that this ratio is one if and only if XD~ = (1 — X)D^, which establishes 
the assertion. □ 
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